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Abstract 

The use of a non-one-dimensional distributed feedback, arising through Bragg diffrac- 



tion in spatially periodic systems (natural and artificial (electromagnetic, photonic) 

H ■ 

crystals), forms the foundation for the development of volume free electron lasers 
(VFELs). The present review addresses the basic principles of VFEL theory and 
describes the promising potential of VFELs as the basis for the development of 
high-power microwave and optical sources. 
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1 Introduction 



Research and development of microwave generators using radiation from an 
electron beam in a periodic slow-wave circuit (traveling-wave tube, backward 
wave oscillator, etc.) has a long history (see e.g. [TTl2f3] ). One of the charac- 
teristic features of such generators is the use of distributed feedback (DFB) 
where the electromagnetic waves produced by the electron beam are collinear 
with the direction of the beam motion (one-dimensional distributed feedback; 
see Fig. [p. 

It is well known that every radiating system is defined by its eigenmodes and 
by the so-called dispersion equation, which describes the relation between the 
frequency to and the wave vector k (u — u(k) or k = k(u)). A thorough 
analysis of the dispersion equation shows that: 

(1) dispersion equation for a FEL p[] in the Compton regime coincides with 
that for a conventional traveling wave tube amplifier (TWTA); 
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Figure 1. 

(2) FEL gain (increment of electron beam instability) in the Compton regime 

1/3 

is proportional to n ' , where n is the electron beam density. 



It was first shown in [5|6] that the law of electron beam instability can change 
significantly under the conditions of a non-one-dimensional distributed feed- 
back formed in a two- or three-dimensional periodic resonator (natural or 
artificial (electromagnetic, photonic) crystal); see Fig. [2j 




diffracted wave 



two-dimensional distributed feedback 




It was shown that the gain of electromagnetic waves changes sharply in the 
vicinity of the points where the roots of the dispersion equation coincide. In 
particular, according to [5|6] , in the case when the DFB in a spatially periodic 
resonator is formed by the two waves participating in the Bragg diffraction, 
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Figure 2. 



the increment of instability is proportional to n^ 4 . If the DFB is formed by N 

number of waves (as a result of Bragg diffraction s = N — 1 number of extra 

waves are produced (see Fig. EJ)), then the increment of instability appears to 

i 

be proportional to Uq +s , provided that the roots of the dispersion equation 
coincide (e.g., for two-wave Bragg diffraction, N = 2 and s = 1, for three- 
wave diffraction, N = 3 and s = 2 ). This result is also valid for electron 
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Figure 3. 

beams moving in vacuum near the surface of a spatially periodic medium, 
in a vacuum slot made inside a periodic medium, or between the diffraction 
gratings of the resonator. The explicit expressions for the starting current j 
were obtained for the conditions of a non-one- dimensional DFB, and it was 
shown that the threshold currents in this case can be sharply reduced. 



The advantages of VFEL are pronounced in a wide spectral range - from 
microwaves to X-rays 



Providing the capabilities for frequency tuning, using wide electron beams 
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(several e-beams), and reducing the threshold current density, required to 
initiate lasing, VFELs offer a more promising potential for the development 
of more compact, high-power and tunable radiation sources than conventional 
electron vacuum devices. 

In particular, experimental investigation of the properties of electromagnetic 
crystals formed by periodically spaced dielectric threads (wires) demonstrates 
that the quality factor of such structures can be as high as ~ 10 6 H-10 8 j9]. First 
experiments have been performed on exciting generation in different types of 
periodic structures (resonators made of diffraction gratings with two different 
periods, resonators based on a photonic crystal formed by metallic wires and 
foils) (for details, see [7]). 

Let us note that in [TT|12|I13] . published some years after [5f6] . the authors, 
considering a particular case of using a two-dimensional DFB for the synchro- 
nization of radiation across a wide sheet electron beam, also came to the con- 
clusion that such distributed feedback can be used for developing microwave 
generators. 

Benefits given by VFELs: 

(1) volume FELs provide frequency tuning by rotation of the diffraction grat- 
ing; 

(2) use of a non-one-dimensional DFB reduces the generation threshold and 
the size of the generation zone. It appeared that if N — 1 + s number of 
waves participate in the formation of a non-one-dimensional DFB, then 
there are N number of connected waves propagating in the electromag- 
netic (photonic) crystal (see FigfS]). This set of N number of connected 
waves has N number of stationary states characterized by wave num- 
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bers ki(u>), (i = 1, ...N). If the relationship (hi — kj)L = 2nmij (jriij are 
the integers, which in the general case are not the same) between these 
waves holds, then a significant change in the threshold current value is 
possible. Let us note that in the case of a one-dimensional DFD (see Fig. 
[1]), these conditions transform into the requirement for the efficiency of a 
one-dimensional DFB [15J; 

(3) VFELs allow the use of a wide electron beam (or several beams) and 
diffraction gratings of large volumes. Two- or three-dimensional diffrac- 
tion gratings (artificial crystals, now often called the electromagnetic or 
photonic crystals) allow distributing the interaction over a large volume 
and overcoming the power restrictions on resonators (see Fig. 2); 

(4) VFELs can simultaneously generate radiation at several frequencies; 

(5) VFELs enable effective mode selection in oversized systems, where the 
radiation wave length is significantly smaller than the resonator dimen- 
sions; 

(6) use of electromagnetic (photonic) crystals with a spatially variable period 
allows increasing the efficiency of lasing [Tfl~4"] . 

2 Radiative instability of beams moving in a spatially periodic 
non-one-dimensional resonator (two- or three-dimensional elec- 
tromagnetic (photonic) crystal) 

Let a relativistic electron beam of velocity u enter the resonator having a 
form of a photonic crystal with length L (the z-axis is perpendicular to the 
crystal surface). Let this photonic crystal have infinite transverse dimensions 
(crystal plate) and lie in the interval (0 < z < L). The set of equations 
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describing the interaction of electromagnetic waves with the "crystal-beam" 
system consists of Maxwell's equations and the equations of particle motion 
in the electromagnetic field. The dielectric susceptibility of a crystal has the 
form e(r; lo) = J2t £ t{u) exp(— iff), where f is the reciprocal lattice vector, e T 
is the Fourier expansion coefficient of e(f, u) and e(f, u) — 1 + w )> where 
x(r, to) is the space-periodic crystal polarizability, x(r, w) = J~J r Xre~ iTr . We 
assume here for simplicity that e(f, u>) is the scalar function. For concretness, 
let us further assume that radiation is excited by quasi- Cherenkov (diffraction) 
radiation (for details, see 05]). We also let \\ T "C 1|. Perturbations of the 
current and charge densities in the linear field approximation may be written 
in the form: 

Sj(k; 00) = exp(— ikf a o) { 5v a (uj — ku) — iu[kdf a (oo — ku)] \ 5n(k; 00) 

a 

= e exp(— ik5f a0 ) ^—i[kf a (u — ku)]} . (1) 

a 

Here 5j(k\ u) and 5n(k; u) are the Fourier transformations of the expressions 
j(f;t) = e^2v a (t)S[f-f a (t)] and n(f;t) = ^5[r - f a (t)) ; 

a a 

u is the unperturbed electron (positron) velocity; 5v a and Sf a are pertur- 
bations of the velocity and radius vectors, respectively, arising due to the 
interaction with the radiation field: 

v a (t) = u + 5v a (t) f a (t) = f a0 + ut + 5f a (t). 
The subscript a denotes the number of the particle. 

Using the equation of particle motion and formula ([1]) for current density, 
one can obtain a set of Maxwell's equations that describes the interaction of 
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electromagnetic waves and particles in crystals: 



k 2 T E{k T i ,oj) — k T '[k T 'E(k T ',u)] E T (k T i,cu)E(k T . 




T 




\7c 2 (u; — k T iu) ) 



f' = 0,fi,f 2 ,..., 



k T = k + t . 



(2) 



The set of equations ([2]) describes a situation when the DFB is formed by 
many diffracted waves (note that in the absence of an electron beam, the 
case of multi-wave diffraction was studied in detail; see e.g. [IE]). However, 
the analysis of such a general situation is very complicated, and so here we 
consider only a two-wave distributed feedback. This allows us to obtain all the 
main characteristics of VFELs analytically and show the advantages of non- 
one-dimensional geometry of distributed feedback over the one-dimensional 
one. 

So let us consider specifically the generation of a cx-polarized wave (the wave 
with the polarization vector orthogonal to the diffraction plane, i.e., the plane 
where vectors k and f lie) for the geometry of the so-called two-beam Bragg 
diffraction [17], where two strong waves are excited, and diffraction occurs by 
the set of crystallographic planes, determined by the reciprocal lattice vector f 
(see FigslHE]). 111 this case, the set of Maxwell's equations, describing two- wave 
diffraction in crystals traversed by the beam, can be written as 




e: = o. 
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(3) 



In equation ([3]), E a — E(k, to) ■ e a , E T a = E(k + r, u) ■ e CT , e CT || [kf], and uj\ = 
47re 2 n /m, where n is the average electron density in the beam. According 
to ([3]), the system "crystal-particle beam" may be considered as an active 
medium with dielectric susceptibility 



• 2 - L (ue^) 2 k 2 c 2 - u 2 



e {k, U) - 1 = £ - 1 o 2 ^7 f^' ^ = = Xt 

7^ 2 7^ 2 c 2 (uo-ku) 2 



col (ue a ) 2 k 2 c 2 — cu 2 
7cu 2 7W 2 c 2 ( w _ £ T #)2 



£o(^r> w ) ~~ 1 — e o — 1 5- 5- 5 — 7^ — -> — ~ 1 e -T — £-t — X-i 



Further, we shall analyze the generation of a wave with wave vector k, making 
a small angle with the particle velocity vector u. In this case, the wave vector 
f is directed at a large angle relative to w, and consequently the 
magnitude of (u) — k T u) cannot become small. As a result, the terms containing 
the expression (u—k T u) in their denominators will be small and can be ignored. 
We shall also neglect the term u\/ r ) - this is justified for real beam densities. 

It is well known that the equation set ((3D has nonzero solutions when its 
determinant equals zero. This defines the dispersion equation, which can be 
written in the form 



(w - ku) 2 (k 2 c 2 - io 2 e )(k 2 c 2 - co 2 e ) - co%e_ T \ = (4) 

„2„2 , ,2w;„2„2 



^(Ue*f( k 2 c 2 _ ^^2 _ u 2 eQ y 



7 cf 



Dispersion equation (111) yields the relationship for uj = u(k) or k = k(cj). 
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3 Generation equations and threshold conditions in the case of 
two-wave Bragg diffraction 



In order to determine the structure of the fields and describe the instability 
evolution in the systems, one needs to know the dispersion equations and their 
solutions, as well as the boundary conditions. 



Now we shall formulate the boundary problem. Let an electron beam with 
mean velocity u be incident onto a plane-parallel spatially periodic plate of 
thickness L. The electron beam is oriented so that it generates radiation under 
Bragg diffraction conditions. Under two-wave diffraction, two fundamentally 
different geometries are possible. 

In the first case (Laue geometry (see FigfJ]), both waves are emitted through 



one and the same boundary of the periodic structure (70 



: kn) 



> and 



7i 



( k T n ) 



> 0, here n is the unit normal vector to the entrance surface). 



Figure 4. Geometry of two-wave Laue diffraction; k, k T are the wave vectors of the 
incident and diffracted waves, respectively, and r is the reciprocal lattice vector of 
the periodic structure. The projections of both wave vectors onto the direction of 
the normal to the surface have the same sign. 
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In the second case (Bragg geometry, see Fig. [5]), the incident and diffracted 
waves leave the plate through the opposite surfaces (70 > and 71 < 0). 























f 






(a 




k 



Figure 5. Geometry of two-wave Bragg diffraction; k, k T are the wave vectors of the 
incident and diffracted waves, respectively and r is the reciprocal lattice vector of 
the periodic structure. The projection of the wave vectors onto the direction of the 
normal to the surface are opposite in sign. 

Let us consider the Bragg case in more detail 

The general solution for the field in a crystal is written as 

4 

E = e a Ciexp(ikir)[l + Siexp(irr)], (5) 

8=1 

where hi is the i-th solution to dispersion equation (jlj) and ki T = ki + T, with 
t being the reciprocal lattice vector corresponding to the planes of diffraction 
reflection. In writing (jSJ), we used four roots of dispersion equation 01]), instead 
of six. Two roots can be discarded because at \\ T \ "C 1, mirror- reflected waves 
can be ignored. 

The boundary conditions, necessary for defining the fields in a crystal, may 
be written as [7115] 

Ci+C 2 + C 3 + C 4 = 1, 
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/ici + f 2 c 2 + /3C3 + /4C4 = 0, 

9\c x + g 2 c 2 + 9303 + g±c± = 0, (6) 
SiCie iXlzL + s 2 c 2 e iK ' 2 * L + s 3 c 3 e iK ^ L + s A c A e iKi * L = 0, 



tu 2 £- T (ue a ) 2 _ k 2 c 2 - co 2 (ue a ) 2 

fc/ T c 2 - (w - kiu) (oj - M) 2 c 2 



The first equation in ([6]) corresponds to the continuity of the incident wave at 
the boundary z = (the wave incident on the crystal is assumed to have a unit 
amplitude). The last equation in fl6]) reflects the fact that behind the crystal, 
no waves move in the diffraction direction (they appear due to diffraction in 
space before the entrance surface of the crystal; the general case, which is true 
for both Bragg and Laue geometries, is considered in [7]). Here i — 1 -r 4 are 
the solutions to dispersion equation @. The second and the third conditions 
correspond to the continuity of the beam density and the beam current density 
at the crystal entrance. Here we apply the expressions obtained from the 
equations of particle motion and the expression for the particle beam current 



ie 2 n (ue a ) 2 k 2 c 2 — u 2 ^ 
mju c 1 (u - ku) 2 

j a = e{ue a )n ^—E a . (7) 

™7 c 2 (u-ku) 



The linear system (Jf)]), defining the coefficients q, has the solution q = A;/ A, 
where A is the determinant of the system and Aj is the z-th minor, obtained 
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clS db result of replacement of the i-th column by 



/ \ 
1 







v / 



Hence, at A — > 0, the field amplitudes inside the crystal increase, and as a 
result, the field in the crystal is nonzero, even when the amplitude of the inci- 
dent wave vanishes. The condition A = with Aj ^ is called the generation 
threshold condition [IB] . Substituting the expressions 



ki = k + kSifi, 
k = u/c, 5i <C 1, 



u — k±u± 



>-0z 



u 7 



(8) 



(where n is the normal to the crystal surface and k$ = (ko z ,k±)) into the 
determinant A, we can represent the generation threshold condition A = as 



(*i - - 5 3 )(5 2 - (5 S ) „ M (6 1 - SJfr - 5 4 )(5 2 - 5 4 ) 



815181 



-s 4 e 



815151 



s 3 e 



(^-^3)(^l-^)(^-^4) 

518181 



s 2 e 



ik5'2L 



(82 - 83) (5 2 - &Q (£3 - £ 4 ) 
515151 



ikSiL 



(9) 







In eq. (Q, the terms containing nonresonant /j and #j (i = 1 -j- 4) were ne- 
glected. 

Upon solving (Q, one may determine the threshold generation conditions, i.e., 
the values of the electron current and other parameters of the beam, at which 
radiation begins to exceed the losses [7118] . 
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In the case of low-gain regime, for instance, the formula for the generation 
threshold under the conditions of two-wave diffraction was obtained in the 
form [5]: 



ir 2 n 2 (ijJi^ 2 



47 



(if) ^ ( X '° + l^l/V^-^" 2 ) (xo + \Xr\/^) sinV x f(y) 



2 



Here a; is the radiation frequency; cu^ = y 4ire 2 no/m e is the Langmuir fre- 
quency of the beam; no is the average electron density in the beam; e and m e 
are the electron charge and mass, respectively; 7 is the beam Lorentz factor; 

= Lu/(un) is the distance traveled by the beam in the crystal; u is the 
unperturbed velocity vector of the beam particles; n is the unit normal vector 
to the crystal surface (directed toward the interior); L is the crystal thickness; 
Xo and Xt are the Fourier expansion coefficients of the crystal dielectric sus- 
ceptibility (their real and imaginary parts are denoted by prime and double 
prime, respectively); (3 = 71/70 is the diffraction asymmetry factor; 70 and 71 
are the cosines of the angles between the normal vector n and the wave vectors 
of the transmitted k and diffracted k + f waves, respectively; ip is the angle 
between the vectors k± and r± (here the subscript _L denotes the projection 
of the vector on the plane perpendicular to u); n is the integer; f(y) is the 
spectral function depending on detuning from the synchronism conditions: 

. (2y + im)siny-y(y + im)cosy 

f(y) = smy ^— , (11) 

y 6 yy + -nn) 6 

where y = kx' 2 L/2 and x-i is the root of the dispersion equation in the absence 
of the electron beam: 



£i, 2 = a±Va 2 + b. (12) 
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Here 



a = -7 (J/70 + k/71) > b = — A (Mt-XtX-t), 

4 47071 

J = O 2 - Xo + T 2 , It = I + a , 

where 6* is the angle between the vectors k and w and a = r(2/c + t)//cq 
characterizes deviation from the exact Bragg condition. 

Let us note that ffTUl) is obtained in the vicinity of the condition imposed on 
the phase difference between the waves in the crystal 

Re (ki z — k 2z )L = 27m, (13) 



which in the considered case has the form 



kWa 2 + b = im. (14) 



and it is assumed that the following inequalities hold: 

% r |L>l, |xr|<l, |xo|<l- 



The condition (ITOl) has a clear physical meaning: the left-hand side of (TlOT) 
contains the term describing generation of radiation by the electron beam, and 
the right-hand side includes the terms describing losses at the boundaries (the 
first term) and absorption losses (the second term) in the medium. It is obvious 
that if the absorption is small, then at fixed values of y (e.g. if y = 0, then 
f(y) = l), eq. (jTOj) yields the following dependence of the threshold generation 
current on the target length (the left-hand side of (TTUT) is proportional to 
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u 2 L ~ n , i.e., ~ j): 



3th ~ (kLj*- (15) 

It also follows from ( flO|) that the value of the threshold current in the case of 
a non-one- dimensional DFB depends appreciably on the parameter /3 and the 
effective photon path length L* in the resonator (recall that L* = Lu/{unj). 



4 VFEL based on an electromagnetic undulator 

Similar considerations enable finding generation thresholds for the beams mov- 
ing in an electromagnetic undulator [T5|2U|21j . 

Let an electron beam move with average velocity uq through a two- or three- 
dimensional periodic medium (electromagnetic (photonic) crystal) placed in 
an external field. Vector potential of the external field can be written as 

(1) A = A w x sin k w z (a magnetostatic wiggler) , 

where A w = B w /k w , k w = 2n/\ w , and x ■ u = 0, with B w being the 
magnetic field strength and X w being the wiggler wavelength (x, y, and 
z are the unit vectors of the Cartesian coordinate system) (c = 1); 

(2) A = A em x sin 

where A em = E em /oj em . Here E em is the amplitude of the electric field 
of the electromagnetic wave and u em , k em are its frequency and wave 
vector, respectively. The external field excites oscillations of electrons, 
having the transverse velocity v± = eA/m / y, where 7 = (1 — v 2 )^ 1 ^ 2 is the 
Lorentz factor of the electron. A moving oscillator emits electromagnetic 
radiation, whose frequency u and wave vector k satisfy the condition of 
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synchronism: 



oj — ku\\ = Q,Q = < 



k w u\\ (for a magnetic wiggler) , 



(16) 



w em — k em u (for an electromagnetic wiggler) , 



where the longitudinal velocity of the electron u\\ = u \1 — J (the 
corresponding longitudinal Lorentz factor I/7 2 = (1 + a w = 

eA w /m is the undulator parameter; k = n(u)u, k em = n(oj em )u em , and 
\n(uu) — 1| <C 1). We shall further consider the case of forward scattering. 



Let us write Maxwell's equation in the Fourier space in the two-wave ap- 
proximation of the dynamical diffraction theory (when the Bragg condition is 
satisfied only for one reciprocal lattice vector). For the field of cr-polarization, 
we have: 



[k 2 - u\l + xo)]E - tu 2 x„ T E T = Aitiuj , 



;i7) 



u 2 Xt E + [k 2 - u 2 {l + Xo )]E T = Amuj T , 



(18) 



where E = E(k, u>) ■ e a and E T = E[k T ,u>) ■ e a are the Fourier components of 
the electromagnetic field of cr-polarization, k T — k + r, e a — k x k T /\k x k T \ is 
the cr-polarization unit vector, j — j(k,uj) • e a and j T = j(k T ,u) ■ e a are the 
projections of the Fourier components of the current density vector onto the 
cr-polarization unit vector, and Xo an d Xr are the Fourier components of the 
periodic tensor of crystal polarizability. 
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Let us use the linear approximation for solving the equations 



In this case, we can write 

f(r,p,t) = f (p) + 5f(f,p,t), 6f~E, <*/</„. (20) 



From ([19]) and flU} follows 
f ^ _)_ # ^ ^ 5 f _)_ ^ 



|+4;W + f ^ = 0. (21) 



In the equation of motion ( 1211) and the expression j = ej vf(f,p } t)dp for 
the current density, we shall go over to the Fourier components. Let us express 
the perturbation of the distribution function in fT2~TT) in terms of E(k,u) and 
then substitute it into (1T7|) and (TTHT) . If the beam's velocity distribution v th 
satisfies the relationships k ■ v t h <C O and k ■ Vth <^ t ■ u, then one can obtain 
the dispersion equation in the form: 



7 

w) = (A; 2 - u 2 (l + X o))(^ - w 2 (l + Xo)) - w 4 *rX-r , (22) 
J=/f> , = (23) 



Let a crystal have a form of a plane-parallel plate of thickness L, then L* 
denotes the particle path length in the crystal (see Sec. Here we shall 
consider only the case of weak amplification, when k" L* <C 1, with k" being 
the imaginary part of the solution to dispersion equation (|22|) . There are two 
characteristic limiting cases for which (1221) can be reduced to the algebraic 
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equation. 

(1) If the inequality 

kvth < 1/L* (24) 

is fulfilled, then we can neglect the beam dispersion and set the initial 
distribution function to be fo(p) = n 5(p — po), where n is the beam's 
spatial density. As a result, we can write the dispersion equation for a 
cold-beam regime. 

D (k,u)(u-kU\\-n) 2 = -2u\k 2 T - uo 2 {\ + xo))F. (25) 

Here F = anc ^ u p = i s the plasma frequency of the beam. 

(2) If inequality ([26]) , inverse to pi]), is fulfilled: 

kv th > 1/L* , (26) 

then we can speak of a "hot" beam regime. To evaluate the integral in 
f[2"5]) . we use the Gaussian distribution function fo(p) = ^°^ e~^~ P0 ^ 2//p *fe 
and obtain the dispersion equation in the "hot" beam regime: 

D (k,u) = -2u\kl - c 2 (l + xo))£- 4= f 1 + *V^^0 , (27) 



fc (for a magnetic wiggler) , 

k ^ em (for an electromagnetic wiggler) 



4-1 Boundary Conditions 



Let a wave vector of the incident wave be k = los. From the continuity of 
the wave field at the boundary follows that inside the crystal, the wave vector 
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— * — * 

k = k + co5n ((5 < 1 because |xo|, \Xr\ ^ 1)- Let us introduce the following 
notations: 7 = k -n/co, 71 = (k +f)-n/uj, & = 70/71, a = ((k +f) 2 -ui 2 )/u 2 , 
Xi — Xo — ct, £ = 7o<5, an d £ = — ^ ~ w anc ^ recas ^ dispersion equations (|2~5|) . 
( |27|) in a more convenient form: 

(e - e?'>)(e - sf')(E - s „) = -(e+ F, (28) 

(e - 4 0) )( £ - 4°') = - (e + lAlf ) 4- ^ f 1 + iv ^^22)) , (29) 



Jo) 



JO) 
£ 2 



4 
1 



Xo - \P1\X1 + ^(Xo+ |/3i|Xi) 2 -4|/3i|XrX-r 



4 L 



Xo - |A|Xi - \/ (Xo + IAIxi) 2 - MPi\XtX-t 



(30) 



In writing (|28[) and (|29[) . we have taken into account that in the case of Bragg 
geometry considered here (the diffracted wave leaves the crystal through the 
entrance surface), /?i < 0. The general solution describing scattering of the 
external wave with the field strength E = E^xe 1 ^^^ by the system " crystal- 
beam" is as follows: 



N 

E = xe i{ ^' wt) Y,( E 3 + e'%)e KjW2/7 ° 



(31) 



and, according to (TX8T) . E T j = E 



-; N = 4 and N = 2 for a "cold" and 



a "hot" beam, respectively. To determine the amplitudes of the fields inside 
the crystal, the boundary conditions must be used. In Bragg geometry, they 
are as follows: 



^1) the continuity of field E at the entrance boundary 



N 



(32) 
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(2) the continuity of field E T at the exit boundary 

SiJTW^ ' (33) 

(3) the continuity of the beam's density at the entrance boundary 
'Sfdp\ s=0 = 0, (34) 

(4) the continuity of the beam's current density at the entrance boundary 
u5fdp\ z=0 = 0. (35) 



Using the expression for the Fourier components df, which follows from (12ip . 
let us express ( 134]) and (135]) in terms of the amplitude Ej in the crystal. It can 
be shown that for the case of a "cold" beam, expressions (154"]) and (135]) yield 
the equations: 

4 E- 4 E- 

These equations, together with equations (132]) and (133]) . give four independent 
conditions that are necessary to define the four amplitudes Ej of the field. 
For a "hot" beam, both of the conditions ( 1341) and ( 1351) can be reduced to 
(|32l . i.e., only the conditions (132]) and (133]) are independent and sufficient 
for defining the two amplitudes of the field. Let g^j denote the elements of 
the fundamental matrix of the system ( 132]) . (133]) . and ( 136]) . The boundary 
conditions enable one to determine the field amplitudes, which can be written 
in the form Ej = E cjji/G, where g~j\ is the cofactor of the element gji, and 
G = det ||<7ii||. If G vanishes, then the amplitudes Ej in the crystal can be 
nonzero, even when the amplitude E of the external incident field equals zero. 
This correlates with the phenomenon of self-excitation (generation), arising 
due to the existence of the distributed feedback. On this account, we shall call 
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the equation G = the generation condition. Calculating the determinant of 
the matrix ||, one can reduce the generation condition to the form: 

J2 \ 1 3 U = (37) 

j= i (ej + j i - .",.) 

for a "cold" beam and 



(3* 



for a "hot" beam. 

^.i? Generation Thresholds 

Recall that the Fourier components of crystal polarizability have the form: 

Xo = -\Xo\ + AXol > X±r = -|X± T | + i\x'± T \ > IXol > IX±rl < IXol > IXrl" 

Now, let us analyze the case when 

V / IXollx / T |«4 0) -4 0) «lx;i- (39) 



As is seen from (13"U|) . this occurs at small deviations from a ~ |^-r(±2y 
|x'ol(l + |/3i|)), i-e., in the vicinity of the absorption edges. The solutions to ( 128)) 
at Re Ei = £ or Re £2 = £, are of particular interest and correspond to the 
synchronism between the beam and one of the DFB resonator modes. In Table 
I, the possible cases are labelled with numbers and letters for convenience. 

In cases la) and 2b), equation (157)) does not have a solution. 
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Table 1 



Synchronism condition 



(0) 



Ree 



Re4 0) = e 



j^(2t/\MM-\x' \(i + \M) 



la) 



lb) 



Viftiix;i-ix(,i(i+ifti) 



2a) 



2b) 



The condition ( l37j) can be satisfied when the following equalities hold: 
in case lb) 



(40) 



in case 2a) 



F(uL,Y 



UJ 



inn 



(41) 



aj 2 a 2 



Let us recall that F = „ 3 P 3 " 2 . The solutions to (!38|) in the case of a "hot'' 



beam have the form 



(1) 



(2) 



u ^g( u W°~ U ^ L * =UJ (\Xo\ ^ 1+ 2^ ~~ \^t\^/\Pi\J l * + 



^n 2 n 2 



n = l, 2... (42) 



a = J-( 2 vfe;|- |xo|(l + |/3i|) 



a 




|/3x| ( 


Uf)L* 


= w ^Ixol" 


n 


= 1, 2... 



2 V / Wix;i-ix , oi(i+iAi) , 



uF ( u (\ "lil+M I "I /hiYW 87rV 



(43) 
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Since F ~ u 2 ~ n ~ j, the value of the threshold current can be found using 
equations (j40|)-fj43l). 

Again, we have jth ~ for a "cold" beam. Equations (|40]) — (T43|) are the 
amplitude conditions of generation. To satisfy the generation conditions, one 
more condition must be fulfilled - the phase condition, common for all the 
cases: 

Re4 0) -4 0) =2im/uL^. (44) 



The condition fj4*4j) leads to the fact that the longitudinal structure of the 
eigenmodes \E\ 2 and \E T \ 2 appears to be close in form to a standing wave. 
Thus, this condition is similar to a well-known condition of standing wave 
formation in a mirror resonator. 

Formulas (|40!) - fj43l) are similar in form. On the left-hand sides, they contain 
the gain of the free electron laser in the appropriate mode. Their right-hand 
sides describe the losses related to energy absorption in crystals and its leakage 
through the boundaries of a DFD resonator. 

It should be noted that in view of the above, the value of the threshold cur- 
rent in the case of a non-one-dimensional DFB depends appreciably on the 
parameter 0i and the effective photon path length in the resonator (recall 
that L* = Lu/{un)). 

Let us consider the equation defining the threshold conditions for backward 
Bragg diffraction (|/3|, = 1 and u || n). In this case, we go over to a one- 
dimensional DFB. The implicit form of the equation for the threshold current 
was obtained in [22] (see eq. (10) in [22]). From this equation, one can also 
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derive the expressions for the threshold current. To do this, let us take eqs. 
(9) and (10), derived in [22], and substitute eq. (9) for the quality factor Q 
into eq. (10). Then we can write 



where n is the integer; V is the dimensionless length of the structure, e.g., 
the waveguide with shallow corrugations (Lf = -L)\ a is the so-called wave 
coupling coefficient due to corrugations, whose physical meaning is similar to 
that of the quantity Xt- Here C = (el n 2 fx/ myco 2 Ng) 1 ^ 3 is the generalized Pierce 
parameter, where / is the current of the electron beam, k is the coupling co- 
efficient between the direct wave and and the electron beam, [i is the electron 
bunching parameter, and N s is the norm of the wave in a smooth wavequide. In 
flUD, $ = -T e L', where T e = -(A + 5). Here A = (u-u B )/u B « [h-h B )/h B 



is detuning from the exact Bragg condition; h B = n/d, u B = cJ g 2 + h%, d is 

the corrugation period, h and g are the longitudinal and transversal wavenum- 
bers, respectively, and 5 is the resonance detuning between the synchronous 
wave and the electrons. 

The function tp'(§) has the form: 



Now let us write the explicit form of the derivative of <p'(&) 

2 2 (-1)" sin $($ 2 - Ti 2 n 2 ) - 4$(1 - (-1)™ cos $) 

At first glance, the spectral functions in ( TTTj) and fH6|) differ. It can be shown, 
however, that virtually they are one and the same function, but of different 



(CL')V($) 



2n 2 n 2 
(aL>r 



(45) 





(46) 
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arguments. Indeed, in the case of a symmetric backward Bragg diffraction, we 
have /3 ~ — 1, and the quantity a in (TT2I) isa = a/4^A. Then assuming that 
5 = and considering the condition (TT4"|) . we get 2y rs $ + to. Substituting 
this equality for the arguments into f|46|) . we obtain ip'(§) = (p'(2y + nri) = 
— ZL T~ f(v)- Thus the condition f|45|) takes the form: 

- — < Ci > "»> = w (48) 

which is equivalent to ( ITUj) (note, however, that in (HSJ), the radiation ab- 
sorption is ignored). Recalling that C 3 ~ /, we immediately obtain that the 
condition ( 1481) yields the same dependence of the threshold current density on 
L as in eq. (flOl . 

We may note in conclusion that in contrast to a one-dimensional case, where 
Laue geometry does not occur, in a non-one-dimensional case it occurs along 
with the generation in Bragg geometry. A more detailed analysis of the radi- 
ation process in Laue geometry can be found in [7] and the references there. 

5 Generation equations and threshold conditions in the geometry 
of three-wave Bragg diffraction 

Above, we have discussed the theory of the volume distributed feedback for- 
mation in the geometry of a non-one- dimensional two-wave Bragg diffraction. 
Let us proceed now to the consideration of multi-wave DFB in a two(three)- 
dimensional periodic resonator. Multi-wave DFB arises when iV number of 
waves simultaneously satisfy the conditions of Bragg diffraction. 

Use of multi-wave Bragg diffraction in a VFEL for the formation of a volume 
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distributed feedback enables one, on the one hand, to appreciably reduce the 
length of the generation area (at a given operating current) or the operating 
current (at a given length of the generation area), and on the other hand, to 
apply electron beams with a large transverse cross section (or several electron 
beams) for generating radiation, which improves the electrical endurance of 
the generator. 

Application of multi-wave diffraction for generating in a microwave range has 
one more remarkable feature - the possibility of selecting modes in oversized 
waveguides and resonators. 

Production of high-power microwave pulses requires high electric strength of 
the generator and radiation resistance of the output window. To reduce the 
load on these elements, the transversal (with respect to the direction of the 
electron beam velocity) dimension of the resonator should be large (much 
larger than the wavelength). As a rule, this leads to a multi-mode generation 
regime and low efficiency When N number of waves are diffracted, the modes 
can be effectively selected due to the requirement to satisfy the Bragg condition 
N — 1 . 

To illustrate the potential of multi-wave distributed feedback, let us consider 
three- wave diffraction in more depth (for details, see Section 13 in [7]). In this 
case, the distributed feedback can be realized in three different geometries: 

(1) Laue-Laue diffraction when the three waves exit through the same surface 
(70,71,72 > 0, here 73 is the cosine of the angle between the wave vector h L 
and the normal to the crystal surface (see e.g. [IS]); Fig. E]); 

(2) diffraction in Bragg-Bragg geometry when 70 > 0, while 71,72 < 0; 
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Figure 6. Three- wave diffraction in Laue-Laue geometry. Projections of wave vectors 
k,k T i, and k T 2 onto the surface normal n have the same sign; t\ and T2 are the 
reciprocal lattice vectors of the periodic structure. 

(3) Bragg-Laue diffraction when 70,71 > 0, and 7 2 < 0; see Figured 




Figure 7. Three-wave diffraction in Laue-Bragg geometry. Projections of wave vec- 
tors k and k T \ onto the surface normal are opposite in sign to the projection of 
vector k T 2\ t\ and T2 are the reciprocal lattice vectors of the periodic structure. 

Similarly to the two-wave case, the problem of the beam interaction with 
the resonator (photonic crystal) can be reduced to the problem of three-wave 
diffraction of an electromagnetic wave incident onto the active medium. The 
active medium here is the system "spatially-periodic structure + electron 
beam" . 

The dependence of the threshold conditions on the length of the interaction 
area at the point of threefold degeneration changes appreciably 




(49) 
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Here in the regime of a "cold" electron beam, G ~ jL 2 , with j being the beam 
current density. Realization of this regime requires the fulfillment of the phase 
condition (k\ z — A^z) L = 2im, (k 2z — k 3z ) L = 27rm, n ^ m. The coefficients 
A and B in (j49l depend on the diffraction geometry, the value of Xt, and 
the indices m and n. According to f|49|) . for a "cold" electron beam, when 
absorption is not important, j thr ~ \MZJ ■ Under dynamical diffraction, 
when the inequality ^ 1 holds, this dependence leads to an appreciable 
reduction in the threshold current. The analysis shows that when a multi-wave 
(iV-wave) DFB is excited, the threshold current depends on L as 



where s = N — 1 is the number of extra waves appearing through diffraction 
in the crystal. 

So, the transition to multi-wave diffraction enables one to significantly re- 
duce the operating current and the longitudinal dimensions of the generating 
system. 

As follows from the above results, the volume distributed feedback (non-one- 
dimensional feedback) has a number of advantages that make its application 
beneficial for generating stimulated radiation in a wide spectral range (with 
wavelengths from microwave and optical to angstrom). Moreover, in a short- 
wave spectral range, where the requirements for the current density and the 
quality of the beam are very strict, it becomes possible to noticeably reduce the 




(50) 



that is, 



1 



Jth 



£3+2s ' 
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threshold current for a given beam propagation area. In this case, the VFEL is 
a unique system providing lasing at relatively small interaction lengths. When 
producing high-power radiation pulses in oversized generators in a microwave 
range, VFELs are beneficial for the reduction of the threshold current, gener- 
ator size, and for mode selection. 

6 Use of a dynamical undulator mechanism to produce short wave- 
length radiation in VFELs 

Numerous applications can benefit from the development of powerful elec- 
tromagnetic generators with frequency tuning in millimeter, sub-millimeter 
and terahertz wavelength range, using low-relativistic electron beams. Low- 
relativistic electron beams in the undulator system can be used for radiating 
in a short- wavelength range, but in this case the undulators period must be 
small. For example, to obtain radiation with the wavelength of 0.3 — 1 mm 
at the beam energy E = 800 KeV— 1 MeV, the undulator period must be 
~ 0.3 — 1 cm. Development of such undulators is a very complicated task. The 
use of a two-stage FEL with a dynamical wiggler generated by the electron 
beam |4J is a possible solution to this problem. 

The above-described potentialities to significantly reduce the threshold cur- 
rents and resonator dimensions by using a non-one-dimensional distributed 
feedback, which arises due to Bragg diffraction in a two- or three-dimensional 
spatially-periodic resonator (electromagnetic (photonic) crystal) and serves 
as a foundation for designing Volume Free Electron Lasers (VFELs), enable 
the development of two-stage VFELs with a dynamical wiggler generated by 
electron beams. A dynamical wiggler can be created with the help of any 
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radiation mechanism: Cherenkov, Smith-Purcell, quasi- Cherenkov [8], or un- 
dulator. VFEL principles offer the advantage of a two-stage generation scheme 
and, in particular, allow one to smoothly tune the period of the dynamical 
wiggler by rotating the diffraction grating. There is a possibility of smooth 
frequency tuning for both the pump and the signal waves by either varying 
the geometric parameters of the volume diffraction grating or by rotating the 
diffraction grating or the beam. Moreover, the VFEL allows one to create a 
dynamical wiggler in a large volume, which is hard to achieve with a static 
wiggler. 

There are two stages in the generation scheme proposed above (see Section 
16, Fig. 16 0): 

(a) creation of a dynamical wiggler in a system with two-dimensional (three- 
dimensional) gratings (in other words, during this stage the electromagnetic 
field, which exists inside the VEFL resonator, is used to create a dynamical 
wiggler). Smoothly varying the orientation of the diffraction grating in the 
VEFL resonator, one can smoothly change the dynamical wiggler parameters; 

(b) radiation is generated by the electron beam interacting with the dynamical 
wiggler, which has been created during the previous stage (stage a). 

Both stages evolve in the same volume. 

7 Conclusion 

Thus, 

(1) volume FELs provide frequency tuning by rotation of the diffraction grat- 
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ing; 

(2) use of a non-one-dimensional DFB reduces the generation threshold and 
the size of the generation zone; 

(3) VFELs allow the use of a wide electron beam (or several beams) and 
diffraction gratings of large volumes. Two- or three-dimensional diffrac- 
tion gratings (artificial crystals, now often called the electromagnetic or 
photonic crystals) allow distributing the interaction over a large volume 
and overcoming the power restrictions on resonators and thus open up the 
possibility of developing powerful generators with wide electron beams (or 
system of beams); 

(4) VFELs enable effective mode selection in oversized systems, where the 
radiation wave length is significantly smaller than the resonator dimen- 
sions; 

(5) principles of VFEL can be used for creation of a dynamical wiggler with 
a variable period in a large volume; 

(6) two-stage scheme of generation, based on a non-one-dimensional dis- 
tributed feedback, can be applied for lasing in the teraHertz frequency 
range using low-relativistic beams; 

(7) two-stage scheme of generation combined with the volume distributed 
feedback can also form the basis for the development of powerful gener- 
ators with wide electron beams (or system of beams); 

(8) use of a hybrid system composed of several phase-locked vircator arrays 
generating modulated beams that enter the VFEL resonator enables fur- 
ther increase of the generator's power output [Tf23] . Moreover, a VFEL 
resonator composed of periodically-arranged metallic elements aids in 
preventing the Coulomb repulsion of electrons. The maximum limiting 
current that can be transmitted through such a resonator is greater than 
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the limiting current that can be transmitted through a vacuum drift space 
of a conventional BWO or TWT (compare with the methods of increasing 
the current transmitted through the resonator, suggested in [2"4"j25f2(j] ). 

To achieve this, two(three)-dimensional electromagnetic crystals formed 
by periodically arranged screens with periodically spaced holes for beam 
transmission can also be used along with those formed by periodically 
arranged wires, wire arrays and cylinders. 

Let us note in conclusion that inverse free electron lasers are known to be 
used for particle acceleration. Similarly, VFELs with a non-one-dimensional 
distributed feedback, arising in electromagnetic (photonic) crystals, can be 
operated in the inverse mode for accelerating particle beams [6]. In this case, 
the acceleration rate is appreciably higher than that achieved in conventional 
FELs. 

References 

[1] J.R.Pierci. Travelling wave tubes (Van Nostrand, Princeton 1950). 

[2] S.H. Gold and G.S. Nusinovich, Rev. Sci. Instrum. 68, 11 (1997) pp. 3945-3974. 

[3] J. Benford, J. A. Swegle, E. Schamiloglu, High Power Microwaves, 2-nd ed. 
(Taylor and Francis Group, New York, London, 2007). 

[4] T.C. Marshall, Free-Electron Lasers (Macmillan Publishing Company, London, 
1985). 

[5] V.G. Baryshevsky, I.D. Feranchuk, Phys. Lett. A 102 (1984) 141. 
[6] V.G. Baryshevsky, Doklady Akad. Nauk SSSR 299 (1988) 19. 



33 



[7] V. G. Baryshevsky, LANL e-print |arXiy: 1101. 0783 [physics. acc-ph, 
physics. optics]. 

[8] V.G. Baryshevsky, K.G. Batrakov, I.Ya. Dubovskaya Journ. Phys. D 24 (1991) 
pp. 1250-1257. 

[9] V.G. Baryshevsky, K.G. Batrakov, I.Ya. Dubovskaya, V.A. Karpovich, 
V.M. Rodionova, Volume quasi-Cherenkov FEL in mm-spectral range, in G. 
Dattoli and A. Renieri (eds.) Free Electron Lasers (Elsiver Science B.V., 1997) 
pp. II-75-II-76. 

[10] V.G. Baryshevsky, K.G. Batrakov, I.Ya. Dubovskaya, Nucl. Instrum. Methods 
A 358 (1995) pp. 493-496. 

[11] N.S. Ginzburg, N.Iu. Peskov, A.S. Sergeev, Pis'ma Zh. Tekh. Fiz. 18, 9 (1992) 
pp. 23-28 (In Russian). 

[12] N.S. Ginzburg, N.Yu. Peskov, A.S. Sergeev, A.V. Arzhannikov, S.L. Sinitsky, 
Nucl. Instrum. Methods A 358, 13 (1995) pp. 189-192 

[13] N.S. Ginzburg, K.E. Dorfman, A.M. Malkin, A.S. Sergeev, Pis'ma Zh. Tekh. 
Fiz. 35, 12 (2009) pp. 9-17 (In Russian). 

[14] V.G. Baryshevsky, A. A. Gurinovich, Experimental study of a volume free 
electron laser with a "grid" resonator, in Proceedings of FEL2006, Aug. 2006 
(Berlin, Germany) p. 335. 

[15] H. Kogelnik and C. Shank, J. Appl. Phys. 43(1972) pp. 2327-2335. 

[16] Chang Shih-Lin, Multiple Diffraction of X-Rays in Crystals (Springer- Verlag 
Berlin Heidelberg New- York Tokyo, 1984). 

[17] Z. G. Pinsker, Dynamical Scattering of X-rays in Crystals (Springer- Verlag 
Berlin, Heidelberg, New York, 1978). 

[18] A. Yariv and P. Yeh, Optical Waves in Crystals (Wiley, 1984). 



34 



[19] V.G. Baryshevsky, I.Ya. Dubovskaya, A.V. Zege, Phys. Lett. A 149, 1, (1990) 
pp. 30-34. 

[20] V.G. Baryshevsky, I.Ya. Dubovskaya, A.V. Zege, Vesti Akad. Nauk. BSSR, ser. 
fiz.-energ. 3 (1990) pp. 49-56. 

[21] V.G. Baryshevsky, I.Ya. Dubovskaya, A.V. Zege, Nucl. lustrum. Methods B 51 
(1990) pp. 368-382. 

[22] V.L. Bratman, N.S. Ginzburg, G.G. Denisov, Sov. Tech. Phys. Lett. 7, 11 (1981) 
pp. 565-567. (Pis'ma Zh. Tekh. Fiziki 7, 21, (1981) pp. 1320-1324). 

[23] V.G. Baryshevsky, A. A. Gurinovich, LANL e-print llarXiv: 0903 .03001 
[physics. acc- ph]. 

[24] D.D. Ryutov, Pis'ma Zh. Tekh. Fiziki 1, 12 (1975) pp. 581-584. 
[25] V.K. Grishin, Zh. Tekh. Fiziki 43, 10 (1973) 2209. 

[26] R.A.Richardson, J. Denavit, M. S. Di Capua, and P. W. Rambo, J. Appl. Phys. 
69, 9 (1991) pp. 6261-6272. 



35 



